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Abstract 

The two-dimensional elastodynamic Green tensor is the primary building block of solutions of linear elasticity 
problems dealing with nonuniformly moving rectilinear line sources, such as dislocations. Elastodynamic 
solutions for these problems involve derivatives of this Green tensor, which stand as hypersingular kernels. 
These objects, well defined as distributions, prove cumbersome to handle in practice. This paper, restricted 
to isotropic media, examines some of their representations in the framework of distribution theory. A 
particularly convenient regularization of the Green tensor is introduced, that amounts to considering line 
sources of finite width. Technically, it is implemented by an analytic continuation of the Green tensor to 
complex times. It is applied to the computation of regularized forms of certain integrals of tensor character 
that involve the gradient of the Green tensor. These integrals are fundamental to the computation of 
the elastodynamic fields in the problem of nonuniformly moving dislocations. The obtained expressions 
indifferently cover cases of subsonic, transonic, or supersonic motion. We observe that for faster-than-wave 
motion, one of the two branches of the Mach cone(s) displayed by the Gartesian components of these tensor 
integrals is extinguished for some particular orientations of source velocity vector. 

Keywords: Analytic continuation. Dislocations, Distribution theory, Elastodynamics, Green tensor. 
Hypersingular kernel. Line source, Mach cone. Regularization. 


1. Introduction 

The two-dimensional elastodynamic Green tensor [1] of the Navier equation is the primary building 
block of solutions for linear elasticity problems involving nonuniformly moving rectilinear line sources, such 
as dislocations (e.g., [2, 3]). Dislocations are the fundamental carriers of plastic deformation in crystalline 
materials [4-6]. Mathematically, a dislocation stands as a discontinuity of the displacement field on its glide 
plane. This discontinuity stands as a boundary condition in traditional methods of solution of continuum 
mechanics [7, 8], or acts more explicitly as a singular source of elastic field if the solution is tackled by 
means of Green functions [9]. In the latter approach, the elastodynamic solution for the strain or stress 
fields involves taking convolution integrals of derivatives of the Green tensor by the source functions [10, 11]. 

Whereas the two-dimensional Green tensor itself is locally integrable, its derivatives are in general hyper¬ 
singular kernels [12], namely, kernels that cannot simply be regularized by means of a Cauchy principal-value 
requirement. Still, they are fully legitimate objects within the theory of distributions, their apparent sin¬ 
gularity being handled by introducing Hadamard’s finite-part prescription [13-18]. From an operational 
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standpoint, finite parts can ultimately be reduced by carrying out suitable integration by parts on convolu¬ 
tion integrals [8]. Hypersingular kernels are commonly encountered in situations involving static [19], as well 
as nonuniformly moving dislocations or cracks [20, 21], notably in the context of so-called boundary-element 
integral approaches [12, 22, 23], and because of their importance in practice, their handling is a recurrent 
issue in wave physics problems [8, 23, 24]. 

In the two-dimensional setting, the dislocation line source is transverse to the (x, y) plane of motion, in 
which it reduces to a point in the idealized case of a Volterra dislocation. In general, point sources lead to 
singular (infinite) fields at the source position and at wave fronts, which poses some problems in numerical 
implementations [8]. However, infinite fields are merely the hallmark of the breakdown of classical linear 
elasticity at the dislocation core. In reality, a physical dislocation has a finite width, that can be measured in 
atomistic simulations or computed by means of specially devised nonlinear models of the cohesive-zone type 
[25-27]. Also, dislocations of a hnite width naturally arise in the framework of gradient elasticity models 
(e.g., [28-32]). 

Being of finite width is a necessary condition for sources to undergo supersonic motion (or faster-than- 
light-speed motion in classical electrodynamics [24]). Indeed, faster-than-wave motions of point sources 
induce Mach or Cerenkov cones with unrealistically infinite field strength [24]. For dislocations or cracks, 
faster-than-wave motion [33] has attracted wide attention during the last decades [27, 34-40]. Also, recent 
medical imaging techniques rely on shear-wave Mach cones induced by a fast moving ultrasonic spot at the 
surface of human skin [41, 42]. Thus, supersonic motion must be allowed for in any comprehensive theory 
of radiation fields generated by moving sources. We should add that, quite generally, the concept of a point 
source can hardly be avoided when no information about the physical nature of the singular source of field 
is available. Then, Hadamard’s finite part regularization, or generalizations thereof, must be employed. We 
refer the interested reader to Ref. [16] for a review of some recent progresses in this direction, motivated by 
the problem of relativistic motion of a point particle in general relativity. 

However, in the specific context of dislocation theory, convoluting the point source by an appropriate 
shape function of finite width that represents the core provides a natural regularization of the relevant field 
integrals at the source location, and at the wave fronts (including Mach cones), and allows one to investigate 
subsonic as well as supersonic motion without the need to address these cases separately [43]. In many 
approaches to finite-size (so-called smeared-out) dislocations [44], core shape functions are often found or 
assumed of power-law decay in the space variable [26, 45, 46]. On the other hand, an exponentially-decaying 
shape function with cut-off characteristic length is produced by the theory of gradient elasticity of the 
Helmholtz type (e.g., [31, 32]) where the convolution is naturally embedded within the Green function of 
the theory as a consequence of the constitutive relations employed. 

This paper introduces an alternative power-law-type way of regularizing the elastodynamic dislocation 
problem, which tames all singularities of the helds in the whole {x, y) plane. While resembling certain 
means [46] currently employed to regularize elastostatic fields in three-dimensional simulations [46, 47] it 
will arise, however, from an immediate analytic continuation of the Green tensor to complex values of the 
time variable, once the elastodynamic fundamental solutions are written down as distributions. Simplicity 
of implementation is indeed a necessary requirement for use in dislocation-dynamics simulations [8, 46]. 

Section 2 reviews several different forms of the two-dimensional elastodynamic Green tensor of the 
Navier equation for the material displacement in an isotropic medium, and its derivatives, which we express 
as distributions. Their regularization is examined in Sec. 3, and applied in Sec. 4 to the computation of 
specific key definite integrals over time, that enter the problem of sources undergoing a velocity jump from 
rest to an arbitrary constant velocity, in the plane-strain and anti-plane-strain settings relevant to screw and 
edge dislocations, respectively. These key integrals —from which expressions of the strain and stress helds 
can be deduced [48]— lead, when employed for faster-than-wave source motion, to Mach cones which we 
further analyze here in terms of distributions. The key integrals are obtained as a difference at their time 
boundaries of non-trivial indehnite integrals. The latter can be used to address more general nonuniform 
source motions since for numerical purposes, a nonuniform motion can in general be represented conveniently 
as a succession of velocity jumps separating small time intervals of uniform motion [27, 49]. The full solution 
for the helds in this problem will be reported elsewhere [48]. A discussion (Sec. 5) closes the paper. 

Notations: Throughout the text, the ‘hat’ notation is employed to denote the unit director a = a/a of 
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a vector a of Euclidean norm a = ||a|| = -^a ■ a. 


2 . Elastodynamic Green tensor and its gradient as distributions 
3.1. Navier equation 

The medium is characterized by the elastic tensor Ciju and the mass density p. The elastodynamic 
Green tensor, Gy, of the anisotropic Navier equation is defined by 

(Sikpdtt - Cijkidjdi)Gkmir,t) = dimS{t)d{r) ( 1 ) 

where d{.) denotes the Dirac-delta distribution and Sij is Kronecker’s symbol. For an isotropic medium, 

Gijkl — T Pi^ik^jl T (2) 

where A and p are the so-called Lame constants. Substituting Eq. (2) in Eq. (1), the isotropic Navier 
equation for the dynamic Green tensor is obtained 

[Sikpdtt -SikP^ - (X + p) didk]Gkjir,t) = Sij S{t)S{r), (3) 

where A denotes the Laplacian. In two-dimensional problems, r = (x, y). The velocities of the transverse 
(shear) and longitudinal waves (sometimes called S- and P-waves) are given in terms of the Lame constants 
as, respectively, 

CT = y/pjp, CL = y/{X + 2p)/p. (4) 


2.2. Anti-plane-strain problem 

In the anti-plane-strain problem, i = j = z. If the material is infinitely extended, the Green tensor 
reduces to Gzz , the usual Green function of the two-dimensional scalar wave equation [50-52], solution of 

{pdtt - pA)Gzz = 5{t)6{r). (5) 


Written as a distribution, its retarded solution reads 


Gtzi^,t) 




( 6 ) 


where 0(t) is the Heaviside unit-step function that restricts the solution to positive times. Here, we have 
denoted with a ‘plus’ superscript the distributional form of the Green function, Gj^, of real arguments, to 
distinguish it from the mere function Gzz(j,t), defined as 




( 7 ) 


and whose arguments will be allowed to take on complex values in the next Section. 
In Eq. (6), the generalized function x^, defined as [13, 14, 18] 


f 0 for X < 0 

( x"^ for X > 0 , 


( 8 ) 


has been used. The derivative of x^^ is given by 

(-f)' = 


( 9 ) 
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_ 2/2 

The derivative of x_f_ ' [13-15] is the following pseudo-function, hereafter denoted by the symbol Pf : 

= ( 10 ) 

Pseudo-functions are particular distributions aimed at regularizing otherwise diverging integrals that involve 
non locally integrable functions (Ref. [13], Vol. I, pp. 38-40; [15], p. 17). They encompass Hadamard’s finite 
part and Cauchy’s principal value prescriptions, and naturally arise upon differentiating certain distributions, 
such as the locally integrable function ' in Eq. (10). 

For notational convenience, we introduce the family of functions 

= (i.<l), (11) 


intended to be continued to complex arguments. They are distinguished from the following corresponding 
distributions of real arguments, denoted with a ‘plus’ superscript: 

/+(r, t) = e{t)9{t - r) Pf(t 2 - = 0 (t) Pf(i 2 _ 

where the symbol Pf is necessary only when v < —2 (note that the positive-time constraint is included in 
this definition). 

Using Eqs. (9) and (10) and the above notations, the gradient of reads 




2/ 2'\~3/2 I Tfc _|_ . . 


(13) 


2.3. Plane-strain problem 

In the plane-strain problem [1, 51], indices in Eq. (3) take on values i,j = 1,2, and the retarded 
distributional solution is the Green tensor {i,j = 1, 2) 


G/Zui) = 


^ f TiTj 


271 p \ 


(t‘ - + (f" - i-voar - c - ’^V4) 


1/2 


,- 1/2 


-{t^-r^fc^r) 


1/2 


(C-rVcaf-C(f^-r74) 


- 1 / 2 ' 


(14) 


It consists of regular distributions made of locally integrable functions; see Refs. [1, 51, 53] for classical 
(non-distributional) writings of this expression. 

The above expression can be re-expressed to emphasize its natural decomposition into spherical and 
two-dimensional traceless (deviatoric) parts. Introducing the traceless tensor 


Ttj{r) = logr = Sij - 2?^ 


t' 


(15) 


the spherical-deviatoric decomposition of the Green tensor reads, in distributional form, {i,j = 1,2) 

\' ^ ;; _L 10„2.,2 .„2^ ^2^0^ 




± //i(u Cpt), 


where a ‘plus’ (resp., ‘minus’) sign applies in the sum when p = T (resp., p = L). The function Gij{r,t) is 
obtained from this expression by simply removing the -I- superscript in /+. 

Using Fourier transforms (FT) proves expedient to derive yet another representation, to be employed 
hereafter. Denoting the FT of G7{r,t) with respect to space variables by Gy (k, t) where k is the Fourier 
wavevector, one has in three dimensions (e.g., [ 21 ]) 


Gij(k,t) 


ejy 

pk 


— sm{cTtk){Sij 
Ct 


kikj) H- sin{ci.tk)kikj 


p 


1 sin(cTtfc) „ 

7 ^ U 

Ct k 


hkj 

p=T,L 


(±) sin(cptA:) 

Cp k^ 


(17) 
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where (±) is a shorthand notation for a factor (±1). Two-dimensional inverses are simply obtained by using 
the convention that k = {k^ + kyY^^ and k^ = k^ = 0. Employing the integrals 



i X cos fp 

2n 


Jo{x), 


^"gcos0e‘"“^^=iJi(x), 


(18) 


where cos0 = k-r, and Jv is the Bessel function, one obtains using Eqs. (6.671-7) and (3.693-1) in Ref. [54] 


d^/c :t.„sin(ctfc) 
(^® k 27r 

A\ 

( 2^)2 


1 1 

— dk Joirk) sin(ctfc) = — (c^t^ — r^) 
iTT Jo 27r 


1/2 




k'^ 


27r 


^ Ji(rfc) sin(ctfc) = — ct 

rC ZTTT’ L ~ 


(19a) 

(19b) 


where in Eq. (19b) the proportionality to Vj of the result was anticipated due of isotropy. Thus, appealing 
to expressions (19a) and (19b) to invert (17) yields in the following alternative form: 


= 

gt (!■> t) = 


^ I ^/^i(cTt, r)5ij + i [d^gj (r, t) + d^gf (r, t)] | 



(bJ = l,2), 


(20a) 

(20b) 


We emphasize that the derivative dig^ is already naturally symmetric with respect to indices i and j, so 
that the explicit symmetrization in Eq. (20a) is merely a matter of convenience for further use in Sec. 4.4. 
Indeed, 


dig+ (r, t) = djgt (r, t) = ^ ^ — 

p=t.l 


TM^(2 cy-r^)- 6 ., 


f\{r,Cpt). 


( 21 ) 


To compute the distributional gradient of C?)*), we observe that = (Jy — rirj)lr. Consequently 
Tij,k{r) = ^ry (r)ffc - ^Tijk{r), (22) 


where we have introduced the totally symmetric and traceless third-rank tensor 


Tijk{r) = djkn + SikTj + SijTk - driTjTk. 


(23) 


According to Eq. (16), is of the type (we omit here for brevity the dependence with respect to time) 
G'+(r) = gs{r)5^j + 5 d(r)Ty (r), (24) 

where gs and gd are scalar functions. Using Eq. (22), the gradient of Eq. (24) reads 


Gy.fe(r) = { g's{r)5,j + 


g'di-f) + -9d{r) 


T’y (r)| ^k - ^gd{r)Tijk{r), 


(25) 


namely. 


Gy-fe(r,i) = ^ ^{/-3(ucpt)[(5y ±ry(r)]rfcT ^(20^12 _r2)/+^(7’,Cpt)ryfc(r)|. (26) 


p=T,L 
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In expanded form, this reads 


G;..(r.*) = + (i‘ - rVcl) 


1/2 

+ 




2 /„2^1/2 
+ 




(e-.V 4 )f-e(e-rV 4 ) 


2 /„2^-l/2 
+ 


nrjTk 




O'Y' C'^jn 




1 


4 ('"-''V 4 ):‘''" + 4 pf(*"-’'“/ 4 ) 


- 1/2 

+ 


-3/2 


L^l 


(27) 


which explicitly features hypersingular kernels interpreted as pseudofunctions, in addition to locally inte- 
grable ones. 

By means of representation (20a) of we can organize terms differently in G^ f.. Using 


dfkdt (r, t) = d%gt (r, t) 


= X! ~ rkflsir, Cpt) - 4(2cpi^ - r^)fliir, Cpt)T,jfc(r) 1 

p=T,L ‘^P ^ J 


(28) 


and 


dkfli{r,CTt) = rkf^3{r,CTt), (29) 

one gets 

+ XI ~ i'fc/-3(u Cpt) - 4(2cp^^ - r'^)f-i{r, Cpt)Tijk{r)\'] 

p=T,L Cp I ^ r j / 

= ^tz,k(^^ 

+ 1^ X ~ - '^d] i’fc/-3(u Cpt) - ^{‘2clt^ - r^)/^i(r, Cpt)ryfe(r) j . (30) 

This expression is employed in Section 4.4. 


3. Regularization 

To motivate the following developments, we point out that in the problem of Volterra dislocations moving 
nonuniformly with time-dependent position s(t) and velocity V(t) = s(t), distributional expressions for the 
material velocity (or the elastic strain field) produced by the dislocations are history-dependent. Namely, 
they typically involve integrals over past times of the following (or of a related) type [48]: 

/ dt'G+^(r-s(t'),<-t')^/(i'): (31) 

J —OO 

where the velocity V(t = —oo) in the remote past is assumed to be a constant (possibly zero), and where 
the upper boundary t' = t is approached by lower values to ensure strictly positive time intervals in the 
derivative of the Green function. The practical necessity for the latter requirement will be illustrated in 
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Sec. 4.2.^ The integral expresses the fact that fields at time t arise as sums of field contributions emitted 
at all ‘retarded’ positions s{t'). An alternative formulation [8] can be derived involving an integral over 
past positions, which requires computing ‘retarded times’. For our purpose, the present setting is more 
straightforward. 

The distributional character of G^- has not been widely acknowledged in the literature, and integrals 
such as (31) have sometimes been dismissed as ill-defined. For instance, it has been advocated that derivatives 
of Gf, be left outside the integral, with expression (31) written in the form 


dk f dt'Gl^{r-s{t'),t-t')Vi{t'). (32) 

J —OO 


Although Eq. (32) is obviously a correct way of carrying out calculations since Gfj is locally integrable, such 
precautions prove unnecessary in the framework of distribution theory where writing (31) is more natural. 
It must be realized, however, that the equivalent expressions (31) or (32) are themselves distributions. 
As will be shown below from Eq. (31), this quantity generates Dirac singularities along Mach cones for 
faster-than-wave motion. 

Although being mathematically well-defined in the sense of distributions and therefore free of non- 
integrable singularities, integrals such as (31) —with arbitrarily prescribed motion s(t )— are inconvenient 
for numerical evaluation. Eor explicit calculations, the Green tensor and its gradient must be regularized. 
The procedure we call hereafter isotropic regularization consists in convoluting G+ by the following isotropic 
representation of the two-dimensional Dirac distribution (Ref. [18], p. 60): 

5(r)=5(a;)^(y) = linr5^(r), with S^r) = (33) 

Here i5'^(r) is a non-singular Dirac-delta sequence with parametric dependence. Eor e finite this corresponds 
to considering a line source with rotationally-invariant core of radius e. We denote this convolution product 

by 

G^ir,t) = [G+ * <5^](r,t) = J dV'G+(r - r',t)S^{r'). (34) 


It should be noted that other regularizing devices have previously been used for dislocations. Eor instance, 
a widely employed model [44, 45] assumes the dislocation to be extended along its direction of motion and 
infinitely thin in the direction normal to its glide plane. In the present context, this would consist in con¬ 
voluting the Green tensor by a one-dimensional Lorentzian core shape —a procedure known as harmonic 
regularization in the mathematical literature [56]— with respect to the coordinate along the prescribed di¬ 
rection of motion. Since it breaks isotropy, this procedure is of considerable complexity in the elastodynamic 
case, for which only the antiplane-strain case (relevant to screw dislocations) has so far been worked out [7]. 
For that reason this type of regularization is not addressed further in the present work. By contrast, the 
regularization considered hereafter allows one to tackle both the antiplane-strain and plane-strain (relevant 
to edge dislocations) cases with only moderate complications. 

Using FTs, convolution reduces to multiplying transforms, and carrying out one Fourier inversion. Inte¬ 
grating over the angle first, the FT of is simply (with Eq. (6.565-3) in Ref. [54]) 


<5^(k) 


I 





Jo{ekz) 

(1 + z2)3/2 


= e 


—ek 


(35) 


^Quite generally, the equal-time value of the time-dependent Green function is defined only as a limit [52] (see p. 189). 
This issue becomes important with the derivatives, and must in general be acknowledged explicitly to carry out calculations 
properly; see, e.g., Sec. 2.2 in [55] for another (related) example where the upper boundary t' = t must be avoided, which is 
achieved there by means of a specific limiting device. Thus, defining the gradient of the Green tensor almost everywhere is not 
enough to get a meaningful result in integrals such as (31). This is another indication —different from the matter of handling 
wavefront singularities at finite time intervals— that the gradient is not integrable in Lebesgue’s sense, but is a distribution. 
We could not find in the literature a proper discussion of this particular point. The present pragmatic treatment, aimed at 
applications, does not pretend to full mathematical rigor. 
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Figure 1: Independent components of in the r = {x,y) plane, at t = 1, for ci^jc^ = 2 (note the changes of scale on 

the G axis). Units are such that cx = 1 and p = 1, and the regularizing width is e = 0.05. The yy component is the same as 
the XX one, rotated by 7r/2 in the (x,y) plane. 


In view of representation (20) of Gj, and since derivatives can be interchanged with convolution, the 
expression of Gf° follows from multiplying in Eqs. (19a) and (19b) the integrands by e“^^, and carrying out 
next the modified integrals. Thus, for instance, 


nOO poo 

/ dA: Jo(A:r) sin = Im / dk jQ{kry 

Jo Jo 


— (e—i ct)k 


= Im[(e — ict)^ + r-^] = Re[(ct + ie)^ — r^] 


where the identity 

\J—z^ = — i sign (Iin z) z for 2 ; S C \ R 

has been used in the last equality. Likewise, the integral in (19b) is modified as 

AU 1 ^_ 

— Ji(fcr-) sin = — Re -[\/(ct + ie)^ — — (ct + ie)]. 

k ^ ' r 



(36) 


(37) 


(38) 


Thus, we arrive at the remarkable result that the regularized form G‘®°(r, f) of the distribution Gfj{r,t) 
is simply expressed in terms of the function Gy(r,t), continued to complex time, as 

Gf° (r, t) = 9(t) Re [Gy (r, t)ct->ct+i e], (39a) 

where our notation means that cxt and cG must be replaced in Gy (r, t) by cxt + ie and cnt + ie, respectively. 
These regularized components are represented in Fig. 1, which emphasizes the smoothness of the wavefronts 
(the value of e employed is arbitrary, chosen for best display). It is noted that G^^ is isotropic in the {x,y) 
plane, equal to twice the spherical part of the plane-strain Green tensor G)|) ; so that, obviously, 

G‘“(r, t) = eit) Re [G,,(r, t),t^,t+u] ■ (39b) 

Analytic continuation has long been known as a method of defining distributions ([14] and [18], p. 159). 
Bremermann’s approach [57-59] consists in defining distributions as the boundary values of analytic functions 
on the real axis, the main variable of integration being extended into a complex quantity. Here, in a multi¬ 
dimensional context involving space and time variables, we make a definite connection between the latter 
approach and the particular analytic continuation introduced here, motivated by the issue of handling a line 
source with finite-size core. 


3.1. Regularized gradient of the Green tensor 

The regularized form of any derivative is obtained from straightforwardly differentiating its regularized 
primitive (e.g., [18] p. 80), which here follows from the commutativity property between convolution and 
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differentiation. Thus, the expression of Gf°i^ results from differentiating Eqs. (39a) and (39b), where the 
functions Gzz and Gy are obtained from Eqs. (7) and (16) by replacing the distribution by the function 

/-i. 

Some care must be exercised when differentiating /_i(r, ct + ie) = [{ct + ie)^ — r^] Appendix A 
emphasizes the fact that whenever the imaginary part of the argument of the square root (computed as a 
principal determination) depends on the variable with respect to which differentiation is carried out, taking 
a distributional derivative generates a Dirac term supported by the branch cut of the square root. In the 
present case however, the group under the square root has imaginary part Im[(ct + ie)^ — r^] = 2ect. As 
this does not depend on r, Eq. (A.lc) tells us that, trivially. 


_d _1__ _Tfe_ 

drk yJ{ct + \£Y -r2 [(ct + i£)2 - r2]3/2 ’ 


(40) 


so that taking spatial derivatives of such algebraic terms is straightforward with the regularization considered. 
It follows that the regularized gradient is simply 


G%°k{r,t) = 6»(t)Re[Gy,fc(r,t)rt^ct+ie], 


(41) 


where Gij^k(j,t) is the function obtained from the distributional expressions of the gradients (13) (anti¬ 
plane-strain) or (26) (plane-strain), by replacing the distributions and by the functions /_i and 

/-3- 



Figure 2: Component near the wavefront cxt^ — = 0, at r = (x, y = 0) and for t = 1, vs. x, for e = 2 with 

fc = 4, 5,..., 15 (solid black). For x < 1, the plots cannot be discriminated when k > 8 (i.e., e < 4.10“®). Units are such that 
Cx = 1 and ^ = 1. Red: representation of the limiting distribution when e —^ 0 (see text). 

Fig. 2 illustrates on Gf°^ the way isotropic regularization handles the finite-part prescription of Pf(c^t^ — 
^2)-3/2 distributional derivatives Gfj ^.(r, t), as e —>• 0. A negative peak of width proportional to e is 

created on the side r > ct oi the singular wavefront = 0. This peak has a sign opposite to that of 

the peak on the other side r < ct. The heights of both peaks are proportional to so that the negative 

peak compensates for the non-integrable part of Gy^fe. In the figure we attempted to represent the limiting 
distribution (27rpcT)~^ Pf(cxt^ — r^)_|_^^^x, symbolizing its finite-part prescription by a downward vertical 
arrow. The imperfect character of this representation illustrates the well-known fact that distributions are 
linear functionals, often unable to deliver useful numbers unless applied to test functions. The only truly 
meaningful plots are those drawn with e finite. 
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4. A case of nonuniformly-moving line source 

4-.1. Fundamental integral, and notations 

We apply the regularized Green functions obtained in the previous Section to the regularization of the 
following fundamental definite integral, defined for t > 0: 

Iijk{r,t) = f G+fc(r - Vt',t- t')dt' = f G+-fc(r -V{t- r),r)dT, 

Jo Jo+ 

{i = j = z, or i,j = 1,2; k = 1,2), (42) 

where the change of variables r = t — t' has been used. This integral arises as a particular instance of 
Eq. (31) with s(t') = Vt', in the problem of a dislocation (or a line force) at rest at negative times and 
instantaneously accelerated to constant velocity V at t = 0 [7, 55]. As the acceleration is instantaneous, this 
is the simplest prototypical case of nonuniform motion. The results obtained below are important elements 
of its full solution [48]. The velocity factor that enters the integrand of integral (31) being a constant in the 
case considered, it has been factored out and omitted in the above definition of Iijk(r,t) (this is the reason 
why, e.g., in Eq. (52) of in Fig. 3 below, a nonvanishing field pattern involving an acceleration wave will be 
obtained for zero velocity). 

For definiteness we focus on the time interval r G (O+jt), but the result will essentially take the form of 
a difference between values of indefinite integrals —denoted by Ik (r, t) (anti-plane-strain case) or Jijk (r, r) 
(plane strain case)— at the boundaries t = t and r = 0+. Thanks to these indefinite integrals, to be derived 
below for any orientation of V, more general time intervals could be considered, which is useful for numerical 
purposes (see Introduction). 

Let c be a placeholder for either ct or cl. The following notations and quantities are employed hereafter: 


f3 = V/c, P=m, 7 = 1/vT^, (43a) 

R(r) =y + \t, (43b) 

S{t) = (43c) 

D(t) = cr —/3 • R(r), (43d) 

H/3) = (1 - - Uifij) ± Uifij, (43e) 


where here and in the rest of the paper the unit vector n = V = /3 indicates the direction of motion, and r 
is a time variable. With respect to this direction, we furthermore introduce the following decomposition of 
the position vector (where rj_ > 0 and r|| can be of any sign): 

r|| = r • n, r'‘- = r-r||n, rj. = ||r*^||, -I-. (44) 


4 . 2 . Anti-plane-strain problem 

In this Section we let c = ct to simplify notations, the longitudinal wave speed cl being irrelevant. 
Going to the co-moving frame by changing r into r -|- Vt, the anti-plane problem consists in obtaining the 
regularized forms of 

Izzk{r + Mt,t)= f G+ ,^(r-f VT,r)dr. (45) 

J 0 + 

With fjL = p<?, the regularized form of this integral is 


I'zzk (r -f Vt, t) = f G'flk (r + Vr, r)dT = Re f -— . U3/2 ^^ 

Jo+ 27rpc Jo+ [[CT -\-ie)^ - 


1 

2TTfl 


■ Re 


^t+ie/c 


/ 0++i e/c 


S{t) '^Rk[T)dT 


r—i e/3 


(46) 
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where the prescription r—J-r — ie/3isa consequence of having shifted the integration path by a value + i e/c. 
It should be noted that the combination of the shifts cr —>■ cr — ie and r —>■ r + ie/3 leaves R = r + Vr 
invariant. Because the imaginary part under the square root S{t) is strictly positive on this path [see remark 
preceding Eq. (40)], the branch cut of the square root is never crossed. Thus, the integral can be computed 
as the difference of boundary values of the following indefinite integral, easily verified by differentiation: 


/fc(r,T) = c 


drS'(T) ^(rfe + Vfcr) 


{rkPp - I3krp){rp + VpT) - CTTk 

(r • A(+) • r) S{t) 


(47) 


It follows that 

Izzki.^ + Vt,t) = Re [/fc(r - ie/3,t + ie/c) - 4(r - ie/3,0+ + ie/c)] . (48) 

This case is sufficiently simple that, using the ‘parallel’ and ‘perpendicular’ notation [Eq. (44)], and 

(r - ie/3) • • (r - ie/3) = [(I - I3'^)r\ + (ry - lePf] , (49) 


we can give the result explicitly as 


/irfe(r + Vt,t) = ^Re 


(1 - /32)r-^ + (ry - lePY 
jrkPp - I3krp){rp + Vpt) - {ct + ie){rk - i£/?fc) 
y/{ct + ie)2 - ||r + Vtp 


irk/3p - /3krp)rp - ig(rfc - ie/Sk) 
i + £2 


(50) 


In the rightmost term, the correct sign of the square root y/ (0+ + i£)2 — = yj—ir'^ + e'^) + i0+ = 

iy/r'^ + £2 stems from the prescription O’*" employed in Eq. (48). This is the reason why approaching the 
upper-time boundary t as a limit, in time integrals such as Eq. (42), is required to obtain a well-defined 
result. This limit must be taken with e finite, before possibly letting £ —>■ 0. 

The quantity 7™fc(r, t) follows from going back to the reference frame by substituting r by r — Vt in this 
expression. Thus, 


Izzkir, t) 2^^ Re ^ - (ct + i £)/3]2 "" 

J {rkPp - PkrpYp - {ct -I- i£)[rfc - {ct + ie)Pk] . {rkPp 
\ y/ {ct -I- i£)2 — r2 


(51) 


Pkrp){rp -ctPp) -i£[rfc - {ct + i£)l3k\ 
yj (r-y - Ctpy + r±^ -I- £2 


Fig. 3 represents t) and '''’Rh £ small enough to highlight the details of the field structure. 

The circular shear wave, centered on the origin, is emitted at t = 0, and expands linearly with time at 
wavespeed ct. The figure illustrates the fact that via analytic continuation, the proposed regularization 
handles faster-than-wave source motion at no additional price. Waves patterns are very different depending 
on the component considered: the source acts as a dipole oriented along Ox in the zzx component, and as 
one oriented along Oy in the zzy component. As a consequence, the Mach cone for |E| > ct in the zzy 
component has branches of opposite signs on either side of the Ox axis, and the leading part of the circular 
wave experiences a change of sign as well. 


4 ..3. Limiting distribution in the anti-plane-strain problem, and Mach-cone analysis 

In this section, we show how to compute the limit of expression (51) as £ —>■ O’*', and investigate some 
of its properties. This limit can be considered a definition of the corresponding distribution 7j^fc(r,t). As a 
warm-up, we consider first the case of zero velocity, for which Eq. (51) reduces to 


^f°fc(i’4)|v=o = 


1 Xk 


— Re 


(ct -I- i £) 


2'Kfl r2 [.^7-2 _j_ £.2 ^_ p2 ^ 


(52) 
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Figure 3: Anti-plane strain problem, with e = 0.01. Components and at time t = 1 vs. position r = {x,y) 

for a line source moving in the positive direction along the Ox axis, at constant velocity V (as indicated) after having been 
instantaneously accelerated from rest at the origin of coordinates at t = 0 (see text for the case V = 0). Units are such that 
cx = 1 and ^ = 1. For better display fields have been thresholded as indicated in the bar legend. 


This expression is nonsingular as r —>■ 0. Obviously, 
1 


lim 


It follows that in the limit, 


= - i(r" - c^t^) 


,2+2^-l/2 

+ 


(53) 


r+ / ,-,1 / 2.2 2\-l/2''’fe 

^)|v=0 — — 2^^ (c i ^2’ 


(54) 


which is locally integrable in two dimensions, and therefore defines a distribution. Indeed, for any function 
/(r) finite at r = 0 and vanishing fast enough at infinity, the following integral (where (p is the polar angle) 
is well-defined: 


/ poo p2T7 

= J J (55) 

Although this is not needed in the above example, we remark that a convenient way of evaluating the 
function 1/r as a distribution (and of obtaining its derivatives) is to write it lim£_>Q 6 {r — e)/r ([18], p. 135). 

We now turn to Eq. (51), assuming that V 0, so that /3 > 0. We must also assume that 1, the 
special difficulty presented by the value /3 = 1 being postponed at the end of this section. To compute the 
limit, we start off by invoking the Sokhotski-Plemelj formula (e.g., Ref. [18], p. 27) 


1 

a; ± i0+ 


1 


= p.v. - 

X 


T i7r(5(x). 


(56) 
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where ‘p.v.’ denotes a principal value. Thus, the overall factor in Eq. (51) becomes 


£->■0 (1 — + [r|| — (ct + ie)/?]^ (1 — + (r|| — ctfiy — 2ie/3(r|| — ct/3) 

= P-V- _ p 2 y 2 ^ I + i7rsign(r|| - ctp)5 ((1 - + (r|| - ctpf) . (57) 

Next, the leftmost square root in the term within braces in Eq. (51) is addressed by means of Eq. (53). 
Besides, because it is of positive argument, the rightmost square root in Eq. (51) is nonsingular; it defines a 
locally integrable term in the limit. Finally, the harmless powers of e in the numerators can be set to zero 
right away. Thus, introducing 

A(r, t) = {y — ctl3) ■ • (r — ctf3) 

= (?'|| - ctpf - {/3^ - l)r5_ (58a) 

= Wr\\ - ctf - (/3^ - l)(r2 - c^t^), (58b) 


one arrives at 


/+ 

zzk 




p.v. 


A(r,t) 


+ i7rsign(r|| — ctf3) 6[A{r,t)) 


X [{rkPp - /3fer-p)rp - ct{rk - ct/3k)] - r^)+^ - 


(59) 

{rk/3p - Pkrp)(rp - ctPp) \ ^ 

V^(r|| - ctpy +r±^ jj 



Figure 4: Locus of the Mach cone for /5 > 1. By definition [Eq. (44)], the Or|| axis is oriented along V and the transverse Or± 
coordinate axis is positive, so that only half of the wave pattern is represented. The Mach cone is tangent at point T to the 
circular wave of radius OT = ct, emitted at t = 0 by the source point M. Given the traveling distance OM = ct/3, the geometric 
construction shows that OP = ct//3. Points on the Mach cone are such that ct//3 < ry < ct(3 and r > ct. Dot-dashed (red): 
supporting line of the Mach cone. 

The function A(r,t) vanishes only at the source position r = Vt = ct/3 if /3 < 1. On the other hand, 
for faster-than-wave motion with /3 > 1 it vanishes for r± = |r|| — ctp}/y/ /3^ — 1, namely, on the supporting 
line of the Mach cone drawn in Fig. 4, which encompasses the source position (because the coordinate 
system employed is such that r_L > 0, only one half of the physical wavefront pattern is represented). Thus, 
Mach-cone contributions for /3 > 1 in Eq. (59) are concentrated in the term ^(A(r, <)). We thus retrieve the 
well-known fact that the Mach cone generated by a point dislocation involves the Dirac distribution [38]. 

We are now in position to simplify Eq. (59). Under the Dirac constraint, expressions (58a) and (58b) are 
employed to express both and in terms of ry. Extracting the real part out of Eq. (59), this 

allows one to turn the terms with square roots in the prefactor of (5(A(r, t)), into rational fractions involving 
absolute values. Factoring signs out of the latter, further simplifications lead to 

Gfe (r, t) = T, 
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where the circular wave C and Mach-cone M vector contributions have, respectively, components 


- 1 f ^ (rk/3p - /3krp)rp - ct{rk - ct^k) 

Ck{r,t)-—{ct -r)^ p.v. , 

Mfc(r,t) = -0(r||,t) (/3^ - 1)+" [(rfe - ctj3k) - PpiPpTk - rp/3k)]S{A{r,t)). 
In definition (61b), the quantity 

^ “ sign(^r|| - ct) sign(r|| - ct/3)] 


(61a) 

(61b) 


(62) 


is a characteristic function equal to 1 in the interval r|| G (ct//3, ct/3), and to 0 otherwise. It thus restricts 
the Mach cone to its physical range (segment TM in Fig. 4). 



Figure 5: Particular velocity directions for which the Mach-cone branch parallel to one of the Cartesian coordinate axes is 
extinguished (dotted) in components Long-dashed, red: directions relevant to component Mx\ short-dashed, green: idem, 

for My. The angle is the half-cone aperture. 


Assuming ,3 > 1, further insight into the Mach-cone part M(r,t) is obtained by using basis orthogonal 
vectors and corresponding coordinates x and y. Let moreover f3 = P(cos(j)ex -b sin (pey), so that 
r||(x, y) = a;cos(^-b y sinc^. The quantity A(r,t) admits the factorization 

A(r,t) = A+(r,t)A_(r,t), (63a) 

(r, t) = (cos 4> ± \//3^ — 1 sin x + (sin cp =F — I cos cpj y — ct(3. (63b) 


Considering t as the main variable in 5(A(r,t)) puts the variables x and y on an equal footing, which is 
most desirable here. Let thus t±{r) be the respective solutions of A±{r,t) = 0. According to a well-known 
property of the Dirac distribution (e.g., [18] p. 49), we have from (63a) 


S{A{r,t)) 


E 

s=±l 


6 {As{r,t)) 
\A-sir,ts{r))\' 


(64) 
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Inserting this expression into (61b) yields 


M(r,t) = — - sign (y cos ^ — xsin^) x 


( 65 ) 


xE 

S=±l 


s COS0 . , \ ^ , Sin 

+ sm 0 Ga; + — COS 0 + s 






'd{r\\,ts{r))S{As{r,t)). 


This expression shows that one of the two branches of the Mach cone component is “extinguished” whenever 
(j) = ±^T or (/) = 7r±^T for M^, and for </) = 7r/2±())T or </) = —7r/2±(/)T for My, where = arctanl/v^/3 — 1 
is the half-opening angle of the cone. For each of these angles, the extinguished branch lies along the Ox 
axis for M^, and along the Oy axis for My (Fig. 5). This accidental extinction, to be further illustrated in 
the next section, is a mere consequence of the tensor character of the integral under consideration. However, 
it is impossible to turn all components off simultaneously, which is physically obvious. 

As to C(r,t), if /3 < 1 the principal-value prescription in Eq. (61a) handles the singularity at the source 
point, which then lies within the circular wave boundary. If on the contrary ,5 > 1, this unique singularity 
splits into two singular points that stand at the intersection of the circular wavefront and of the Mach cone 
lines (tangency points represented by T in Fig. 4). Their investigation can be carried out straightforwardly 
by decomposing the overall rational fraction into partial ones, using the factorization formula (63a). As 
these field singularities are well known [38, 39], the analysis will not be pursued further. 

Before closing the section, we return for completeness to the marginal case /? = 1. Then, Eq. (57) is 
meaningless, because we face the problem of a double root in the denominator. The case must be addressed 
by the following generalization ([14], p. 60) of the Sokhotski-Plemelj formula: 


1 

(a:±i0+)2 


= Pf ^ ± i 7rd'(x) 


( 66 ) 


which stems from differentiating (56). Thus for /3 = 1 equation (57) becomes 


lim -— 

E->o [r|| 


1 

{ct -I- ie)]^ 


= Pf 


1 

(r|| - ctf 


i7r^'(r|| — ct) , 


(67) 


The analysis can then be completed as above, separating t) into “circular wave” and “cone” parts. 

By (67), the cone part is formally proportional to (5'(r|| — ct) = —c~^S'{t — r^^/c). It vanishes, as must be, 
due to the identity ([18], p. 36): 


f(t)S'{t - to) = f(to)5'{t - to) - f'{to)S{t - to), 


( 68 ) 


where / is some differentiable function at t = to- Applying this identity with to = ''’II/c to the prefactor 
of d'(t — r||/c) —considered as a function of time— shows that the whole term cancels out. The surviving 
circular-wave contribution reads 




2^-I 

+ 


Pf- 


{rk/3p - l3krp)rp - ct{rk - ctjSk) 


(r|| - ct)2 


(/3 = 1 )- 


(69) 


We note that this expression is nonsingular even when r_L = 0. Indeed, the numerator of the fraction reduces 
then to 


{rkl3p - Pkrp)rp - ct{rk - ct/3k) = (ry 


ct) Tk - (rii -I- ct)Pk 


{I3 = l,rx = 0). 


(70) 


Moreover, [c^t^ — = {c^t^ — ry)^A^ go that the distribution behaves as |ct — r||| -which remains 

regularized by the Pf prescription. 
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4-.4- Plane-strain problem 

Evaluating the plane-strain integral 

rt 


Iijk{r -\-Vt,t) = [ G'^+ (r-1-VT,T)dr, i,j,k=l,2, 

Jo+ 


(71) 


is notably more complicated. We employ the same method as for the anti-plane-strain case, our efforts being 
directed towards obtaining the key indefinite integral that allows for a direct solution. 

The form (20a) of and the associated representation (30) of its gradient carry over mutatis mutandis 
to a similar representation for In terms of the latter the integral we need to compute is: 

liTki'^ + Vt, t)= [ G)"°fe(r -H Vt, r)dr 


/o+ 


= i::,ir + Yt,t)S, 


+ ‘ Re V W r"" i 

•7o+-Hie/cp 2 




S{t) 


R{rY 


+ {i^ j) 


(72) 


r->r-i£/3(P> 


where = Y/cp, and where 7‘™j,(r-f Vt, t) is given by Eq. (48). By the same arguments as in the previous 
Section, the above integral follows from the indefinite integral 


J^Jkir,T) = c j Y 


^ik 


Sir 


Rj (' t ) 


+ ^ j) 


= cUr ^ 


R(rY 

ry(%))-<5, 


2c^r^ — RirY 


2SiTf L ^ ‘"J RirfSir) 

in which the wave speed c is generic. Let us introduce the notation 

'^iiki^^^) ~ Jijk{r,t)c^cp 


Tijk (R-('J')) 


(73a) 

(73b) 

(74) 


to denote the tensor Jijkir, r) in which c is substituted by the specific wave speed Cp with p = T, or p = L. 
Expressed in terms of j)j))(r, t), Eq. (72) reads 

ifYkir+ Yt,t) = P::^{v+ Yt,t)5,, 


+ ^ 4^fc(i’-ie/3^^\i + ie/cp)-4'’fe(r-ie/3^P\o+-bie/cp) 

^^P p=T,L S ^ 


(±) r ap). 


(75) 


Directly carrying out integral (73b) to obtain Jijk is of considerable difficulty. In view of Eq. (73a), we 
instead compute it as the second symmetrized gradient of the following simpler indefinite integral that we 
verify in Appendix B: 


Jj (r,T) = c j dr Sir) 


Rj jr) 

RirY 


= ^ + 7 log[ 7 -D(r) -b S'(T)]r™A^^^(/3) - r^L/(T)B^ij(n)| 

where is defined in Eq. (43e), Bmiji^) is the third-rank tensor 
Rmtj i^) — -f i^ziripP^pP^l A ^zlpP^p'am) dj2j 


TIT, 


(76a) 

(76b) 

(76c) 
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and Li{T) is the vector 


L(r) 


/ 

V 


log- 


— arctan 


R{r) 

r ■ R(r) 

r X (3)^S{t) 


\ 

) 


(76d) 


In the latter expressions, is the Levi-Civita symbol, and (r x = ezki^kPi- Equation (76b) holds 
up to an arbitrary time-independent and r-dependent integration ‘constant’ in the right-hand side, -which 
explains why the latter blows up as /3 —0 (i.e., V — 0 ). We call such arbitrary quantities time-independent 
terms (TITs). Contributions from TITs are ignorable, since they ultimately cancel out in the difference 
of boundary values in (75). For the sake of completeness, we provide in Appendix C another form of 
(not used hereafter), in which the integration constant is adjusted to yield the expected finite limit 
as V — 7^ 0. 

The contraction LiBuj in (76b) can be made more transparent by computing the components of Buj. 
One finds that Bmin) = —r;i(ri), B 2 i 2 {n) = 7]i(n) and Bii 2 {n) = i? 2 ;i(n) = —r; 2 (n). Thus, introduc¬ 
ing the following traceless matrices that correspond to pure-shear and simple-shear modes of plane-strain 
deformation (e.g., [60]): 


; -1 ) - ( 

0 1 \ 

. 1 0 j’ 

(77) 

one has 



i?.i,(fi) = - E = 

i?.2,(fi) = - T2^{n)ulf = 

a^l,2 

^zipTpj (n), 

(78a) 

(78b) 

and 



LiBuj=- LiTiUn)U^f. 


(79) 


a=l,2 


We can now obtain J^-fc as the second symmetrized gradient of Jj{r,T). In view of definition (73a) of 
Jijk by an integral, this way of proceeding involves interchanging the integral and the space derivatives in 
this definition. Strictly speaking however, 

^ [dfkJji^, T-) + d%J,{Y, r)] = Jyfc(r, r) -f Sing., (80) 

where we emphasize the fact that the left-hand side is equal to Jij-fc(r, r) only up to time-dependent singular 
Dirac-like distributional parts concentrated on the source path. Indeed, Eq. (D.4) in Appendix D shows that 
such a Dirac term, of support the line = 0, arises when differentiating the arctangent in (76d), because the 
latter function is a discontinuous function of r. This can be traced to the fact that the integrand in definition 
(76b) of Jj is inversely proportional to the norm of R(t) = r -|- Vr, which vanishes at time t = —r • V /V 
when = 0 . Since R is invariant under the complex-valued shifts considered in our calculations [see 
remark following (46)], this kind of singularity is not removed by analytic continuation. Thus, for instance, 
the integral 

/*t+i e/c 

c dr 

JO+S-ie/c 


S{t) 


Rji-r) 


R{t)\ 


(81) 


whose second derivative enters (72), is ill-defined (infinite) on the segment C{t) = {r j r = —tV,t G (0,t)} 
contained in the glide plane r**- = 0, which represents the trajectory of the source expressed in the co-moving 
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frame. Quite generally, the source path is a notorious source of difficulties (see, e.g., [8] and references 
therein). 

Here however, the difficulty is somewhat artificial, as it stems from the route adopted to reach Jyfc(r, r). 
For brevity, we bypass the problem by remarking that the latter is necessarily a function without any 
Dirac term, since by (73b) its time-derivative is the integrand, which is free of Dirac-like singularities. It 
follows that extracting Jijk(r,T) from a double differentiation of Jj{r,T) can be done simply by dropping 
any Dirac-like distributional parts that would normally arise in the process of differentiating the arctangent 
in Eq. (76d). 

Then, completing the calculation reduces to carrying out cumbersome differentiations. Introducing new 
vectors 


Q = R/5, q = r/(r.A(+)-r), (82) 

one has 

CT 

diS = -Qi, dij \og{jD -\-S) = {qiPi - l3iqi)Qi “ + TIT, (83) 

where an ignorable term in the second expression has been lumped in the TIT notation. In terms of these 
gradients, the first symmetrized derivative of Ji reads 


I 


-b jj,f) = + PjQ^) + \A[J{qjPi - fdjqi) + A^Jiq.Pi - fiqi)] r^Qi 


CT 

^ L 


Aim + Aj^ qi rm 

+ 27log(7D -|- S)A^j ^ — 2 LmBimj — rm {Li^iBmij + + TIT, 

where gradients of the components of vector L defined by (76d) are as follows: 


(84) 




CT Ri 


ffiBy — 


^zip 


CT Rp 


(/3-Q-f)gp 


(85a) 

(85b) 


In arriving at (85b) the following identity has been employed: 
[(r X l3)zS]'^ -b (r • R)^ = (r • • r) R^. 


( 86 ) 


Unless the TIT is suitably chosen in Eq. (76b) (which we shall not bother to do), is not a symmetric 
tensor. This lack of symmetry is our motivation for having symmetrized Eq. (20a) in the first place. 
Introducing further 


Uij — dij -b QiQj^ 


in terms of which 


u, = 


A, 


(-) 




r ■ A(+) • r ’ 

Wk^ 


VF,, = - 2 


r 

'’m.-n-mi '3 
r • H(+) • r 


dkQi — „Uki, dkPi — aC-i-I 

b r • A*^+J • r 


(87) 


( 88 ) 
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the second symmetrized derivative (80), amputated from its singular part, finally provides 


Jijk{T^,T) — 2^2'I g{UkjUkif3' 


;) 


where 


+ [{Wk^Pl - mkl)Vmj + {WkjPl - ^JWkl)V^^] rmQl 
+ [{nPi - l3iri)Vkj + {rjPi - Pjri)Vki + 2{rkPi - l3kri)Vij] Qi 

+ ^ [Vim(rj/3l - Pjri) + VjmiTiPl - j5iri)\ TmUkl 


CT 


g “1” {Yik^^j “t“ VjkTi -j- ^VijTk') 


{Li iBkij ~t" LijBkii + 2Li kBiij') {Li ik^^mij ~t“ Li jkBjnii^ ? + TIT, 


Lx,ik — 


CT 


Ly ik — 


(Uik — 2RiRk^ , 

Lx,pk + (Ukl/3l 'gQk'^ Qp + (j3 ■ Q ~ 


CT\ Wkp 


A(+l • 


+ TIT. 


(89) 


(90a) 

(90b) 


The zero-velocity form of Jijk{r,T) is most easily obtained, not from Eq. (89), but rather by differentiating 
twice expression (C.l) of Jj(r, t)|v=o- One finds 


</yfc(r,T)|v =0 



TiTjCk - ^\/c^ 


- r2 ryfc(r) -f TIT. 


(91) 


This expression blows up when r —^ oo, but this is not a problem. Indeed the part that explodes, namely 
— (c^T^/r^)Tyfc(r), cancels out in expressions such as (75) where it is eliminated [after division by the 
relevant wave speed squared; see (75)] from adding the longitudinal and transverse contributions, which 
have opposite signs. 

Substituting Jyfc(r, r) and -I- Vt, t) as given by (48) into Eq. (75) completes our solution for 

/‘®^(r -I- Yt,t). The quantity follows from going back to the reference frame by means of the 

substitution r —r — Vt. The analysis of the distributional limit £ —>■ 0 is too long to be reported, but would 
follow the same lines as in the antiplane-strain case. 

Fig. 6 illustrates these plane-strain calculations on the component for supersonic motion at 

velocity |E| > cl, with focus on the branch-extinction property. The outer circular wave is of longitudinal 
character; the inner one is transverse. Various branches of either Mach cone (or both) can be extinguished, 
or almost extinguished, on any component of the tensor Rjk for specific orientations of the velocity. In the 
figure the fundamental critical angles for Mach-cone branch extinction are (^t = arctan(/3^ — 1)“^/^ ~ 0.411 
(shear cone) and = arctan(,5£ — 1)“^^^ ~ 0.923 (longitudinal cone) radians. Eor an angle such as </> = tt/S 
close —but not very close— to these critical angles or to one of their symmetric counterparts (see previous 
section), a simultaneous quasi-extinction of branches on both cones is observed. 


5. Summary and discussion 

To summarize, we emphasized the often overlooked distributional character of the Green tensors of the 
plane-strain and anti-plane strain components of the Navier equation, and their gradients, and studied a 
regularization procedure for these quantities. The regularized expressions are suitable to direct exploitation 
in practical problems involving nonuniformly moving source lines. The regularization consists in spreading 
the point source over a finite width £, by means of an analytic continuation of the elastodynamic Green 
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Figure 6: Plane strain problem. Component for supersonic velocity V = 2.5ct, and various orientations of the 

velocity vector (as indicated). Units are such that cx = 1 and /r = 1, cl = 2cx, and s = 0.01. 


tensor to complex time. The connection between this procedure and one of convolution by one particular 
source shape function has been established. 

As an application, some important definite integrals of tensor fields over finite time intervals were obtained 
from non-trivial indefinite integrals, for both plane-strain and anti-plane-strain conhgurations. These definite 
integrals, to be used elsewhere [48], are encountered in the problem of a moving line source suddenly 
accelerated from rest to a constant velocity. In this context, they determine the material velocity and 
strain fields. Numerical computations were carried out to illustrate the fact that, since the Green function 
embodies all wave propagation effects, cases of subsonic, transonic, or supersonic source motion can be 
handled altogether by one single closed-form expression. Moreover, a property (of geometric origin) of 
Mach-cone branch extinction on Cartesian components of the tensor integrals under consideration was 
reported. 

There are two ways of envisioning the obtained integrals. The first one is to consider them as computa¬ 
tional intermediates to reach the fields in the form of distributions. Indeed, the method is straightforward 
since we do not need to care about causality constraints at wavefronts. As was explicitly demonstrated 
in the anti-plane case, the distributional forms are obtained by letting the source width parameter e —>■ 0 
(Volterra limit). In this limit, the resulting distributions represent helds generated by the moving point 
source. However (see Introduction), such distributional Volterra solutions are devoid of any physical mean¬ 
ing in the vicinity of their singularities, namely, the source position and the wavefronts. This is particularly 
obvious for faster-than-wave motions where Mach cones show up in the limit in terms of the Dirac distri¬ 
bution. Thus, no meaningful numbers can be extracted from the limiting integrals, unless the latter are 
convoluted with a suitable source shape function to produce physical fields. Of course, using for this purpose 
the isotropic shape function employed for regularizing the Green tensor from the outset would result in the 
same expressions with e finite reported above (assuming the calculation to be manageable this way, which 
is far from granted; this is why we started from Green’s function instead). Other shape functions could in 
principle be employed, perhaps at the price of additional difficulties. 

The alternative point of view —which has our preference— is to consider that the isotropic shape function 
we employed is of physical relevance (at least approximately, say, as some effective mean core shape; e.g., in 
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the sense of an angular averaging), and to use for simulations the obtained definite integrals as they stand. 
This spreads Mach cones over a width of order e, and makes fields finite everywhere. Keeping e nonzero 
reveals that dramatic variations in strength and sign of physical fields take place near wavefronts, as our 
figures clearly show, which may be of consequence for the dislocation nucleation problem [8]. 

Then, there remains the question of the value to be attributed to e. From a purely numerical standpoint, 
any value not leading to overflows in field expressions is acceptable. On the other hand, physics requires 
both the source position and its width to be prescribed by coupled equations of motion, which makes e{t) 
a function of time [27]. For dislocations in metals, this “physical” e does not exceed a few multiples of the 
interatomic distance (but can also be notably smaller, depending on the velocity and acceleration regimes). 
From this perspective, the definite integrals obtained above must be though of as a post-processing device to 
estimate fields in the surrounding medium, once the dislocation position and “physical” width are known as 
functions of time; the typical size e considered in the present work being substituted by the “physical” one.^ 
However, it is important to point out that our present derivations do not account for the specific radiative 
contributions induced by such a time-dependence of e. With applications in mind, it could nonetheless be 
feasible to tentatively neglect the latter (assuming it to be subdominant) and carry out this substitution 
within each time interval in a numerical procedure of the type alluded to in the Introduction. 
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Appendix A. Diflferentiation in the complex plane in the sense of distributions 


We argue here that differentiation in the sense of distributions of principal determinations of logarithms 
or square roots in the complex plane must be performed bearing in mind that the derivative involves a distri¬ 
butional part. The latter is supported by the negative real half-axis, on which the principal determinations 
undergo a jump in the direction parallel to the imaginary axis. Thus, 


dlog 2 =-h 2\'k9[— Re2;)5(Im2;)d(Imz), 

z 

A-sfz = + 2 i(— Re 2;)y^(5(Im 2 ;)d(Im z), 

2 y z 

dz-i/2 = -iz-3/2d2- 2i(-Rez);^/^(5(Imz)d(Imz). 


(A.la) 
(A.lb) 
(A.lc) 


As a rule, differentiation, with respect to some variable m, of the logarithm or square root of a complex¬ 
valued function z = f{u) creates a Dirac-like distributional part, if: (i) the imaginary part of f{u) depends 
on u, and (ii) there exists values of u such that Re/(it) < 0 while Im/(it) = 0. For instance, from Eq. 
(A.lc) follows the non-trivial result that 


d 1 

dt ^(^ct -I- ie)2 — 


c{ct + ie) 

[{ct + ie)2 — 7'2j3/2 



S{t) 


(A.2) 


(this equation is not needed in the main text). 


^Ideally, the regularizing source shape should, for consistency, be the same as that employed in deriving the equations of 
motion of the source. This would be asking for too much in the present state of knowledge (the equations of motion of Ref. 
[27] were obtained within harmonic regularization). 
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Appendix B. Proof of Equation (T6b) 


The proof proceeds by differentiation. One finds, in succession 

Id D 
= s- 

c dr 

-;^log[ 7 ^(r) +5'(r)] = 1/(75'), 
c dr 

Id CT — S{t) 1 r • R 
cd7 R{t) 

1 d 

- -r-a] 
c dr 

It follows that 


— Jfc(r,T) = 


r • R(r) cr 

(r X I3)z 

(r X (3)zS{t) S 

i?2 ■ 

m n4;\(3) 


(3‘^S (3‘^S 

/32S'i?2 i 

DPk nA\-\f3) 

a2Q ' 02 Q 

P2QU2 


[(r • R)Tifc(n) - cr 

^zrnn'I^rriPn^zip'^pkij^^] i 


where the last form results from Eqs. (78). Moreover, since indices m^n^i^p take on values 1,2, 


^zmn^zip — ^im^np ^mp^in- 


Therefore, 


l3^SR^-^Jk{Y,T) = \Dj3k + riA[fJ{l3) \ - n [(r • R)rife(n) - ct{rij3p - rpj3i)Tpk{n)]. 
dr L J 

A straightforward expansion of the right-hand side then shows that 

P'^SR^^Ju{v,T)=P^S^Rk. 


whence the result. 


(B.la) 

(B.lb) 

(B.lc) 

(B.ld) 


(B.le) 


(B.2) 

(B.3) 

(B.4) 

(B.5) 


Appendix C. Alternative indefinite integral for the plane-strain case 


By going to the limit V —>■ 0 in the integrand, the integral in (76a) naturally reduces to 
Jj (r,T)|v=o = c J dr a/ c^r^ - 


r 

^2 L' 


cr\/ c^r^ — 


log ^cr + \/c^r^ — 


Const. 


(C.l) 


The arbitrary implicit time-independent constant in (76b) can be adjusted to retrieve this limit, by employing 
the same equation, but now with 

1 CT - S 1 ,0X2 

log— - -logr- -(r •/3) 


/ 


L = 


R 


— arctan 


r • R 


(r X f3),S 


■ sign(r X /J)^ 


(C.2) 


It should be noted that as the regularization considered in the paper involves adding to r an imaginary 
part proportional to (3 [see Eq. (75)], the quantity (r x (3)z is always a real number. Although it is more 
regular (the expression of Ly being now continuous as /3 —>■ 0), the above expression of L differs from (76d) 
by ignorable TITs. 
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Appendix D. Distributional derivative of the arctangent in L 


We clarify hereafter the issue of the distributional derivative of the arctangent in L. First, it is recalled 
that arctana; is a continuous function, while axcian{l/x) is discontinuous at a; = 0. The discontinuity 
manifests itself in the following well-known identity: 


1 TT 

arctanx -I- arctan — = — sign a:. 

X 2 


It follows that, in the distributional sense. 


(D.l) 



TT 

— Sign X — arctan x 



1 

1 -I- a:^ 


(D.2) 


This distributional derivative differs from the usual one by the presence of the Dirac term in the right-hand 
side, and plays an important role when working with single-valued material displacements generated by 
dislocations [21]. 

Consider now h{x) = arctan]/(a;)(a;)] where both functions / and g change signs at their zeros. In 
view of the above derivative, a little reflection shows that the derivative h'(x) must be computed in the 
sense of distributions as 


. d fix) 

h {x) = — arctan 

da; g{x) 


_d_ 

da; 


^ sign]/(a;)] sign]g(a;)] - arctan 
2 J[x) 


f'[x)g{x) -g'{x)f{x) 
p{x)+g‘^{x) 


Trsign[f{x)]g'{x)6{g{x)), 


(D.3) 


where the subscripts indicate that differentiations are carried out with g and / kept constant, respectively. 
The Dirac peaks in the derivative at the zeros of g{x) are observable numerically when h'{x) is approximated 
by [h{x + g) — h{x — g)]/{2g), with g a small number. 

Applying this formula to the differentiation of Ly in (76d) yields 


= -^ {{qx I3)^{2S + r ■ Q)Qi - [cr(q x /3)^/3i -f q • Re^i^/3m]} 

— TT sign(r • R ){ez^m|3m)5{{r X f3)z). (D.4) 

Up to TITs, the regular part can be reduced to (85b), while the Dirac part, supported by the line r-*- = 0, 
is not time-independent since it involves R = R(t) in its prefactor. This is the source of the difficulty 
encountered in Section 4.4. 
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